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Введение. Статья посвящена решению задачи оптимального 
управления процессом получения биогаза при непрерывном 
режиме его переработки в двух метантенках. Цели работы: 
представить математическую модель данного процесса и 
разработать алгоритм выбора оптимального управления. 
Материалы и методы. Созданная математическая модель 
описывает получение биогаза из отходов животноводства 
при последовательной переработке субстрата в двух метан- 
тенках. Рассматриваются случаи одинаковых и различных 
температурных сред (мезофильной и термофильной). Для 
данной модели сформулирована задача оптимального 
управления в виде задачи Лагранжа. Ее управляющими 
параметрами являются скорости поступления субстрата в 
метантенки. Предложен алгоритм решения данной задачи, 
основанный на чис 

ленной реализации принципа максимума Понтрягина. При 
оптимизации применялся гибридный генетический алго- 
ритм с дополнительным поиском в окрестности лучшего 
решения методом сопряженных градиентов. 

Результаты исследования. Разработана новая математиче- 
ская модель, описывающая процесс получения биогаза при 
последовательной переработке субстрата в двух метантен- 
ках. Предложен и программно реализован численный алго- 
ритм решения задачи оптимального управления.Численные 
исследования показали, что для термофильной среды ско- 
рость образования биогаза практически вдвое выше, чем 
для мезофильной. Установлено, что последовательная пе- 
реработка субстрата в двух метантенках с одинаковыми 
температурными средами позволяет вдвое увеличить ско- 
рость образования биогаза. Если температурные среды в 
метантенках различны, то в первом из них следует исполь- 
зовать мезофильную среду, а во втором — термофильную. 
При этом скорость образования биогаза несколько ниже по 
сравнению со случаем, когда в каждом из метантенков ме- 
зофильная среда, однако степень переработки субстрата 
выше на 10—15 %. 
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Обсуждение и заключения. Полученные результаты могут 
быть использованы для расчета и конструирования биога- 
зовых установок, а также при разработке соответствующего 
программного обеспечения. 
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Введение. При работе животноводческих предприятий, в частности, птицефабрик, свинокомплексов, 
ферм крупного рогатого скота (КРС), образуется большое количество отходов. В результате их сбраживания 
(ферментации) в специальных резервуарах (метантенках) можно получить топливный газ (биогаз) и ценные 
органические удобрения [1]. Этот процесс называется метаногенезом. Он может протекать при периодическом 
и непрерывном поступлении субстрата. В первом случае метантенк наполняется однократно и полностью опо- 
рожняется по завершении ферментации. Во втором случае одновременно и непрерывно происходят два процес- 
са: поступление субстрата и удаление его переработанной порции. 

Обычно для метаногенеза используется мезофильная среда (при температуре 25-38 °С) или термо- 
фильная (45—60 °С). Оптимальная температура для мезофильной среды — 37 °С, для термофильной — 56 °С. 
Продолжительность ферментации для этих сред равна 25 и 12 суткам соответственно. 

Экономическая эффективность производства биогаза зависит от различных факторов: вид и количество 
сырья, климатические условия [2] и т. д. Кроме того, на выработку биогаза существенно влияет скорость по- 
ступления субстрата в метантенк. Значение данного параметра зависит от объема метантенка и вида сырья. 
Чтобы найти оптимальное значение указанной величины, требуется решить задачу оптимального управления. 

Следует отметить, что при непрерывном режиме ферментации субстрат не успевает пройти полную пе- 
реработку. Для повышения выработки биогаза необходимо использовать два метантенка, чтобы субстрат по- 
следовательно перерабатывался в каждом из них. Различные аспекты такого процесса изучены в работах [3, 4], 
а его технические реализации отражены в патентах [5-7]. 

Во работах, посвященных математическому моделированию метаногенеза, например [8, 9], не пред- 
ставлен поиск оптимального управления либо находится его асимптотическое значение [10-12]. Аналитическое 
решение задачи оптимального управления получено в [13]. Однако описанная там математическая модель су- 
щественно отличается от представленной в данной статье. Предлагаемый здесь численный метод решения за- 
дачи оптимального управления применим к широкому классу моделей. 

Материалы и методы. Схема последовательной переработки субстрата представлена на рис. 1. 


Метантенк 1 Метантенк 2 
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Рис. 1. Схема двухстадийного режима метаногенеза 





Пусть [0 — концентрация питательных веществ в субстрате, поступающем в первый метантенк 
(кг/м°); [ и Г› — концентрация питательных веществ в субстрате для первого и второго метантенков соответ- 
ственно (кг/м’); Хр и Х› — концентрация метанообразующих бактерий в первом и втором метантен- 
ках (кг/м°); О; — объем субстрата в первом метантенке (м), О› — во втором; 40 — скорость поступления 
субстрата в первый метантенк (м”/сут), 41 = "2 и 42 = скорости выделения субстрата из первого и 


второго метантенков соответственно. 
Потерь субстрата не происходит, поэтому 
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40 = 41 = 92 =Ч4. 


Величина Чо определяется по объему отходов, количеству и производительности биогазовых устано- 








ВОК. 
Относительные скорости поступления субстрата для метантенков составляют: 
та 1 а 
те и ‚902 _ 92 
Оо @ О Ш 0 
Объемы субстратов в метантенках равны 
Ч 4 _ щ О, 
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Для описания метаногенеза при последовательной переработке субстрата используется математическая 
модель, основанная на математической модели популяционной динамики метанообразующих бактерий для од- 
ностадийного режима метаногенеза [11, 14]: 
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Здесь нижние индексы переменных ( Х,, [,) и параметров (Нис г, Нид г, а, В, Вг, и; ) соответствуют 
номеру метантенка Ге {1, 2}. Параметрами модели служат Ии.; И Иша; — максимально возможные относи- 
тельные скорости соответственно прироста и отмирания бактерий (сут '); В; — безразмерный коэффициент 
усвоения субстрата; а; и Ь; — эмпирические коэффициенты (м\/кг); ^ — параметр, равный нулю, если темпе- 


ратурные среды в метантенках различны, и единице, если эти среды совпадают. 
Значения параметров модели определяются в соответствии с выбранным температурным режимом ме- 
таногенеза (мезофильным ИЛИ термофильным). Если температурные среды в метантенках совпадают, то 


Итет = Ите2 = Ито › Ита1 = Ита2 = Ита ИТ. Д. 
Начальные условия: 
Х! (0) = Х, (0) = Хо, 11 (0) = 15 (0) = 1%, 
где Хо — естественная концентрация метанообразующих бактерий в исходном сырье; Гу равно концентрации 


питательных веществ в непереработанном субстрате. 

Модель (1) построена в предположении, что в метантенках поддерживается оптимальная и неизменная 
температура процесса. Если температурные среды различны, то при поступлении субстрата из одного метан- 
тенка в другой должен происходить нагрев или охлаждение до температуры во втором метантенке. 

Скорость образования биогаза (м”/сут) в -м метантенке равна 


УИ те ИХ 
и; =—— , 
а; +2; 
где у; == коэффициент, характеризующий скорость преобразования питательных веществ субстрата В биогаз 


(м. м/кг). 
Для получения задачи оптимального управления необходимо дополнить систему уравнений (1) крите- 
риальным функционалом 





Т 
&Х [5 Х. 
= | УИ 1 72Нте? 2—2 @-> тах, (2) 
0 


“1 +Ё 92 +2 
определяющим общий выход биогаза с 1 м* субстрата в первом и втором метантенках за время Г. Оптимизи- 
руемыми параметрами задачи служат относительные скорости поступления субстрата в метантенки и и и В 


системе уравнений (1). 
Системы (1)-(2) являются задачами Лагранжа. В общем случае задачи оптимального управления име- 
ется система дифференциальных уравнений вида 


(хе), (0), 0), го, т жек вая”, (3) 
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ее начальные условия: 


х (0) =хХо. 
Требуется найти оптимальное управление и (г), доставляющее максимум критериальному функциона- 
лу 
Т 
7=[ю(<(9, ч(1), 6) тах. (4) 
0 
Для численного решения задачи на промежутке [0,7] вводится разностная сетка с узлами 1 =0, 4, 
2, В, И, -. и =Т с постоянным шагом й=1;.1 —& [15]. 


Численное решение системы дифференциальных уравнений (3) осуществляется методом Рунге — Кутты 
четвертого порядка: 





Хх =Х; + +2К> +2Кз +К4), 1=0, 4—1 › 
К = м, в), 


й й 
к [ме мьи+ (5) 


й й 
К: =Ё х,+—.К, м, Е +- |, 
3 [ 2 2 "] 


Для решения задачи удобно ввести 4^-мерный полный вектор управлений Ч = (и;), 1=а4. При этом 


разностной аппроксимацией критериального функционала (4) является выражение: 
И 
1=/(0)=У 8 (хр, ц,,)-й > шах. (6) 
1=0 


Задача решается в приведенной ниже последовательности. 

1) Задается полный вектор управлений Ц. 

2) Численно решается исходная система дифференциальных уравнений (3) с помощью соотношений (5) 
и вычисляется значение критериального функционала (4) с помощью разностной аппроксимации (6). 

3) Численно интегрируется система сопряженных уравнений (в направлении «справа налево») согласно 
соотношениям: 





ОР : 
рр =рни+й- +0; рн" 8=Ъа-1, (7) 


1 


ра =0. 





В ОЕ (х,, и, й,) 
Здесь р = р(#) — двойственные переменные принципа максимума Понтрягина; С; = Эх матри- 
1 
ца Якоби, составленная для системы (3). 

4) Выполняется процесс оптимизации по вектору 1. В данной работе в качестве метода оптимизации 
использовался генетический алгоритм с вещественным кодированием и дополнительным поиском в окрестно- 
сти лучшего решения с помощью метода сопряженных градиентов. 

Для рассматриваемой задачи оптимизации метаногенеза вектор фазовых переменных равен 


х = соот (ХТ, [, Х>, 12) ‚ вектор управления и = союп (м, и> ) ‚ а матрица Якоби имеет вид: 




















Ни Ники я ата ; Ант 0 0 
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Результаты исследования. Фиксируются теоретические значения переменных, найденных в результа- 
те численного решения системы дифференциальных уравнений математической модели метаногенеза. Учиты- 
ваются отклонения указанных значений от экспериментальных значений относительно вектора параметров [10, 
16] (табл. 1). Значения параметров модели метаногенеза оцениваются путем минимизации суммы квадратов 
этих отклонений. 

Таблица 1 


Оценки значений параметров модели метаногенеза для мезофильной / термофильной сред 






































Источник сырья 
Коэффициент 

Птицефабрики Свинокомплексы Фермы КРС 

Ито 0,482 / 0,821 0,346 / 0,783 0,297 / 0,563 

Ита 0,353 / 0,528 0,291 / 0,423 0,254 / 0,351 

а 34,781 / 43,875 7,242 / 21,653 5,013 / 8,733 

Ь 116,457 / 14,674 37,347 / 9,278 18,722 / 5,455 

В 2,344 / 3,189 1,495 / 2,084 1,413 / 1,983 

у 1,463 / 1,963 1,373 /1,907 1,299 / 1,813 








Рассматриваются следующие способы переработки субстрата. 
1. Используется единственный метантенк. 


П. Используются два метантенка с последовательной переработкой субстрата, температурные среды в 


которых совпадают. 


При этом для каждого метантенка возможны два варианта среды: мезофильная ИЛИ термофильная 



































(табл. 2). 
Оптимальные параметры метаногенеза 
ое ристиаа Источник сырья 
Птицефабрики Свинокомплексы Фермы КРС 
Т. Один метантенк с мезофильной / термофильной средой 
и”, сут" 0,149 / 0,402 0,128 / 0,341 0,112 /0,268 
и’, м’/сут 7,25 / 24,51 3,58 / 10,76 1,81 /5,42 
В: 56 / 43 51/38 52/41 
П. Два метантенка с мезофильной / термофильной средой в каждом 

и, сут ' 0,168 / 0,432 0,144 / 0,362 0,117/0,287 
и> ‚сут 0,491 / 2,448 0,412 / 1,858 0,367 / 1,424 
и, м’/сут 15,19 / 50,21 7,79 /22,19 3,85 / 10,92 
12/1, % 45 /35 38 / 33 38 / 33 














Ш. Два метантенка, в первом — мезофильная, во втором термофильная среда / 


в первом — термофильная, во втором — мезофильная среда 




















и] ‚сут ' 0,149 / 0,402 0,128 / 0,341 0,112 / 0,268 
из, сут" 0,312 / 0,051 0,222 / 0,059 0,183 / 0,051 
и, мУсут 15.35 /25,12 6,61 / 11,42 3,26 / 5,42 
2 / Ц, % 30/34 29/35 26/32 
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* * 
Здесь и —_ оптимальные значения относительной скорости поступления субстрата; иж — соответ- 


* 


ствующий ИМ среднесуточный выход биогаза; —_ —_ полнота переработки питательных веществ субстрата. 
0 


Графики скорости образования биогаза показаны на рис. 2. 


Скорость образования биогаза: 
— — в схеме с одним метантенком 
в схеме с двумя метантенками 











Скорость образования биогаза: 
— —_ в схеме с одним метантенком 
в схеме с двумя метантенками 








0 20 40 60 80 100 
+ сут 
6) 
Рис. 2. Скорость образования биогаза при использовании отходов свинокомплексов для двух метантенков 
с мезофильной (а) и термофильной (6) средой в каждом 


Ш. Используются два метантенка с последовательной переработкой субстрата, температурные среды в 
которых отличаются. 

При этом рассматриваются два варианта. 

1) В первом метантенке мезофильная среда, а во втором — термофильная. 

2) В первом метантенке термофильная среда, а во втором — мезофильная. 

Оптимальные значения параметров см. в табл. 2. 

На рис. 3 приведен график скорости образования биогаза. 
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у, Скорость образования биогаза: 
— — в схеме с одним метантенком 
м` (сут в схеме с двумя метантенками 

















0 20 40 60 80 100 
Ь сут 
а) 
и’ Скорость образования биогаза: 
— — в схеме с одним метантенком 
м` /сут в схеме с двумя метантенками 


Винни 

ИГ. 
ДТ 
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1, сут 











6) 
Рис. 3. Скорость образования биогаза при использовании отходов свинокомплексов и двух метантенков: в первом 
метантенке мезофильная, во втором — термофильная среда (а); в первом метантенке термофильная, 
во втором — мезофильная среда (6) 


Использование одного метантенка при непрерывном режиме ферментации позволяет переработать пи- 
тательные вещества субстрата на 40-45 % для мезофильной среды и на 50—55 % — для термофильной. 

При использовании двух метантенков с одинаковыми температурными средами доля переработки суб- 
страта увеличивается на 10 % для мезофильной и на 5 % для термофильной сред. При этом среднесуточная вы- 
работка биогаза возрастает практически вдвое. Скорость образования биогаза для термофильной среды практи- 
чески вдвое выше, чем для мезофильной. 

Два метантенка, эксплуатируемые в различных температурных средах, демонстрируют различные ре- 
зультаты. Так, при использовании первого варианта доля переработки субстрата примерно на 25 % выше, чем 
при использовании одного метантенка с мезофильной средой. При втором варианте доля переработки субстрата 
на 5—10 % выше, чем при использовании одного метантенка с термофильной средой. Кроме того, организация 
процесса по первому варианту обеспечивает практически двукратное увеличение скорости образования биогаза 
(если сравнивать с одним метантенком с мезофильной средой). Во втором варианте значимый прирост скорости 
не отмечается (если сравнивать с термофильной средой). Из вышесказанного следует, что использование второ- 
го варианта нецелесообразно. 

Первый вариант несколько проигрывает в скорости образования биогаза случаю использования двух 
метантенков с мезофильной средой в каждом, хотя обеспечивает более высокую (на 10-15 %) степень перера- 
ботки субстрата. Это связано с тем, что при использовании двух метантенков с мезофильной средой оптималь- 
ные значения относительных (а следовательно, и абсолютных) скоростей поступления субстрата выше, т. е. в 
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единицу времени перерабатывается больше субстрата. Скорость образования биогаза для второго варианта 
вдвое ниже, чем для случая использования двух метантенков с термофильной средой. 

Обсуждение и заключения. Разработаны математические модели, описывающие метаногенез при по- 
следовательной переработке субстрата в двух метантенках, когда температурные среды в них совпадают и от- 
личаются. Этим моделям соответствует задача оптимального управления процессом метаногенеза на основе 
принципа максимума Понтрягина. В работе представлен алгоритм ее решения. Управляющими параметрами 
служат относительные скорости поступления субстрата в метантенки. 

В результате численных исследований обнаружено, что последовательное использование двух метан- 
тенков с одинаковыми температурными средами позволяет увеличить степень переработки субстрата на 5-— 
10 %. При этом скорость образования биогаза возрастает вдвое. Степень переработки субстрата для мезофиль- 
ного режима также на 5—10 % выше, чем для термофильного. При этом скорость образования биогаза в термо- 
фильной среде практически вдвое выше, чем в мезофильной. Это связано с более высокой интенсивностью 
процесса (выше оптимальное значение относительной скорости поступления субстрата). 

Также показаны результаты эксплуатации метантенков с различными температурными средами. В пер- 
вом варианте субстрат сначала перерабатывается в мезофильной среде, а затем поступает в метантенк с термо- 
фильной средой. Во втором варианте, наоборот, субстрат из метантенка с термофильной средой поступает в 
метантенк с мезофильной средой. При первом варианте перерабатывается больше субстрата. Если сравнивать с 
применением одного метантенка с мезофильной средой, преимущество составляет около 25 %. Если сравнивать 
с использованием одного метантенка с термофильной средой — 5-10 %. Для первого варианта скорость обра- 
зования биогаза практически вдвое выше по сравнению со случаем одного метантенка с мезофильной средой. 
Во втором варианте скорость образования биогаза почти такая же, как при использовании одного метантенка с 
термофильной средой. Таким образом, применение второго варианта нецелесообразно. 

Для первого варианта скорость образования биогаза несколько ниже, чем в случае применения двух 
метантенков с мезофильной средой в каждом, но в то же время степень переработки субстрата выше на 10- 
15 %. Скорость образования биогаза во втором варианте вдвое ниже, чем при использовании двух метантенков 
с термофильной средой в каждом, что еще раз показывает неэффективность второго варианта. 
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